The very long range nature of capillary interactions in liquid films 
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Micron-sized objects confined in thin liquid films interact through forces mediated by the deformed 
liquid-air interface. This capillary interactions provide a powerful driving mechanism for the self- 
assembly of ordered structures such as photonic materials or protein crystals. Direct probing of 
capillary interactions requires a controlled force field to independently manipulate small objects while 
avoiding any physical contact with the interface. We demonstrate how optical micro-manipulation 
allows the direct measurement of capillary interactions between two micron sized spheres in a free 
standing liquid film. The force falls off as an inverse power law in particles separation. We derive and 
validate an explicit expression for this exponent whose magnitude is mainly governed by particles 
size. For micron-sized objects we found an exponent close to, but smaller than one, making capillary 
interactions a unique example of strong and very long ranged forces in the mesoscopic world. 



It is a well known fact that small objects floating on a 
liquid surface cluster together. Bubbles on the surface of 
a soap solution Q or cereals in a bowl of milkQ attract 
each other with long range forces arising from the in- 
terface deformation under particles weight (or buoyancy 
for bubbles). Close packed configurations for such macro- 
scopic objects are found to minimize gravitational poten- 
tial energy. Shrinking lengths to the mesoscopic scale, 
particles weight soon become too weak to produce any 
significant deformation and hence attraction. However if 
the particles are confined in thin liquid films a deforma- 
tion of the interface is unavoidable. This is the case, for 
example, of a colloidal suspension drying on a solid sub- 
strate, or dispersed in a free standing thin film. When the 
thickness of the liquid film becomes smaller than the bead 
diameter, the interface has to deform with an increase 
in surface energy. The liquid interface will then react 
on the particles with forces aiming to reach a minimum 
surface (energy) configuration, that usually corresponds 
to close packed two dimensional crystals |3j- Such phe- 
nomena, already observed by Perrin in 1909 [4], have at- 
tracted considerable interest in recent times due to their 
relevance for the engineering of photonic materials Q and 
protein crystallography [6]. Consequently, a strong effort 
has been devoted to the theoretical analysis of the in- 
volved forces, resulting in a long series of papers reviewed 
inQ. Prediction for macroscopic objects have been con- 
firmed by experiments on immersed cylinders [8[ or par- 
ticles attached to holders 0. 

However no experiment so far has been able to directly 
measure the strong capillary force acting between an iso- 
lated pair of mesoscopic objects, despite the fact that 
it is in the mesoscopic and nanoscopic realm that this 
effect finds the most interesting applications. Any phys- 
ical contact with the particles would inevitably produce 
a significant deformation of the liquid-air interface and 
dramatically affect the interaction. On the other hand, 
due to long ranged hydrodynamic interactions in 2D, par- 
ticles mobility are very sensitive to interparticle distances 



FIG. 1: (Color online) Wetting geometry. Complete wetting 
is assumed, so that the liquid film wraps the particle inside a 
spherical cap of radius r c . For r > r c the liquid- air interface 
is free-standing and slowly falls to the large distance height h. 
The local height of the interface measured from h is indicated 
by the function C( r )> whose gradients are assumed to be small 
everywhere (ip c <C 1). 



and force measurements are difficult to deduce from par- 
ticles trajectories. A static, highly non-invasive method 
is required for a direct and reliable measurement of these 
interactions. 

In this Letter we demonstrate how optical micro- 
manipulation [lQ] allows the precise measurement of cap- 
illary interactions between two micron sized spheres con- 
fined in a free standing thin liquid film. Holographic 
optical tweezers [ll| allow to isolate a single pair of par- 
ticles and scan their relative distance from close contact 
to several tens of diameters. Capillary forces will tend 
to push the particles towards each other and out of the 
optical traps until the restoring trap forces will balance 
the attractive force. The intensity of capillary force can 
be then deduced by measuring particles displacements 
from trap centers, after trap calibration. Working in a 
free standing liquid film is essential for accurate capillary 
force determinations since no particle-substrate interac- 
tions have to be taken into account. 

The force law can be anticipated by calculating the sur- 
face tension forces acting on a pair of spherical objects 
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confined in a thin liquid film. Kralchewsky et al.[l2[ de- 
rived the shape of the meniscus around the two particles 
using the method of matched asymptotic expansions. An 
implicit expression for the force law was obtained, whose 
evaluation required the numerical solution of a system 
of nonlinear equations in the accessible system param- 
eters. Confining here to mesoscopic objects, we show 
that Kralchewsky solution for the force can be very well 
approximated by an inverse power law whose parame- 
ters have an explicit expression in terms of the system's 
physical properties. We assume complete wetting (zero 
contact angle) with the liquid surface departing from a 
contact ring of radius r c with continuous slope angle i(j c 
(Fig. la). Surface tension will exert a force on the con- 
tact ring whose resultant is orthogonal to the ring plane 
and has a modulus: 



F = 7 2i\r c sin ip c 



(i) 



where 7 is the liquid- air surface tension (7 ~ 35 mN/m 
in our experiment). For small gradients r c = asin^ c ~ 
a^c, z c = a(l — cos^ c ) ~ a?/r?/2, and the force F ~ AtvjZc 
is then proportional to z c with a strength of order 1 
nN/nm. The order of magnitude of other forces into 
play is 1 pN for particle's weight and 100 pN for the 
maximum optical force exerted by our trap. The inter- 
face is supposed to be flat and horizontal far from the 
particle and we indicate with £(r) the local vertical dis- 
placement of the interface from this reference surface. In 
the small gradient approximation, Laplace equation for 
the pressure drop across free standing portion of the up- 
per surface (r > r c ) reads [l3|: 



7 V 2 C(r) = pg((r) 



(2) 



where p is the liquid density and g the acceleration 
due to gravity. The only axisymmetrical solution to ([2|) 
vanishing at infinity takes the form [3]: 



bottom surface will have a symmetrical shape to the top 
one. Therefore two equally strong capillary forces will 
act on the top and bottom contact rings of an isolated 
particles but they will cancel out giving no net force. 
When a second particle is introduced in the film at a 
distance s, it will in turn contibute to the interfaces de- 
formation producing a tilt of the contact lines around 
the first particle (Fig. lb). For small deformations we 
can still think of contact lines as circles but this time 
slightly tilted. Within the superposition approximation 
by Nicolson[14j, the amount of tilt would simply be given 
by the gradient of the surface deformation produced by 
an isolated particle at the location of the second one. 
According to (|4j) the interface height around an isolated 
particle is expected to decay logarithmically with the dis- 
tance producing a tilt of the contact lines that falls off 
as the inverse interparticle distance. Capillary forces F 
acting on the contact rings will not be balanced anymore 
but a net attractive force would appear whose intensity is 
given by the projections of F on the film reference plane: 



f(s) = 2F(s) 



d((r) 



dr 



16^72^(5) : 



(5) 



where s is the interparticle distance. The resulting 
force is than proportional to the square of z c and in- 
versely proportional to the distance separating the parti- 
cles. In analogy with Coulomb electrostatic force in 2D, 
the quantity z r is usually called the "capillary charge" of 
the particle[16j. However, the capillary charge is not an 
intrinsic constant property but a slowly varying function 
of particle separation. This function can be calculated 
by imposing the continuity and differentiability of the 
interface across the contact ring. 



C(r c ) + C(s) +h = a- 
r c = a sin ip c ~ \J1az c 



(6) 
(7) 



C(r) 



t&nxpc K (qr) 
q Kx(qr c ) 



(3) 



where Ki{x) is the modified Bessel function of ith 
order [15] and q~ l = ^j/pg is the capillary length. This 
is the length scale below which gravitational forces play 
no role in determining the interface shape and it's about 2 
mm in typical solvents. Therefore we can safely replace 
Kq(x) and K±(x) with their small argument expansion 
and write for small gradients: 



C(r) = -2z c [log(gr/2) + 7e 



(4) 



where j e is the Euler-Mascheroni constant. At the mi- 
cron scale gravitational energy will be negligible and the 



where we rely on the superposition approximation [14J 
to express the interface vertical displacement field as the 
sum of two single particle displacement fields (|4j). We 
are also assuming in (|6]) that the deformation field £(s) 
produced by one particle is constant over the contact line 
of the other. Equations I6|7I can be solved analytically 



giving: 



z c {s) 



h 



W [-q 4 s 2 a(a - h) exp[4 7e - l]/8] 



(8) 



Where W is the Lambert-W function [17[. The above 
expression for z c is a slowly varying function of s in meso- 
scopic systems (a < s <C q~ l ) and it can be well approx- 
imated by its logarithmic expansion about s = 2a: 
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z c (s)~z (s/2a) a 

z z c (2a) 
d\og z c 



dlogs 



2z 



s—2a 



h 



(9) 
(10) 

(11) 



Accounting for changes in z c the attractive force ([5|) 
will still display a power law behavior but with an expo- 
nent smaller than one: 



-2a 



(12) 



where a is an explicit function of the three lengths: 
particle size, film thickness, capillary length of the sol- 
vent. The exponent 1 — 2a is equal to 1 for h — > a but 
then quickly drops to a fairly constant value as soon as 
h is small enough to produce significative forces (~ 1 
pN ~ 10 z KbT/o). This constant value depends practi- 
cally only on the particle size and varies very little when 
using different typical solvent properties. Changing par- 
ticle size from 10 nm to 10 jiva produces a corresponding 
exponent variation in the range 0.92 to 0.82. For our 
particle size we predict an exponent of 0.86. 
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FIG. 2: (Color online) Schematic view of the experimental 
setup described in text. 

A schematic view of the experimental setup is reported 
in Fig. 2a. 

Latex beads (5.7 jim diameter) are dispersed in a 2/3 
water - 1/3 glycerol mixture with added 0.2% wt surfac- 
tant (SDS). A thin film is obtained by sweeping the so- 
lution on a square frame (10 mm side) of nylon wires (60 
fim thickness). The free standing liquid film is enclosed 
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FIG. 3: (Color online) Film thickness profile around an 
isolated particle, b) Monochromatic light reflected off the 
film displays interference fringes centered around an isolated 
trapped particle. A period in the fringe pattern corresponds 
to 240 nm thickness variation. Figure a) reports as open cir- 
cles the extracted film height f (r) versus distance from parti- 
cle center r together with a logarithmic fit (red line). 



in a humidity chamber and placed over the 40x objec- 
tive of an inverted optical microscope (Nikon TE2000- 
U). The same objective is used to focus the laser beam 
(A=532 nm) reflected off a spatial light modulator (Holo- 
eye LCR-2500) into two, dynamically reconfigurable, op- 
tical traps [ll|. Axial confinement is guaranteed by the 
normal components of capillary forces. We can also ac- 
cess the thickness variations around an isolated particle 
by viewing the film under reflected monochromatic light. 
To this end a red diode laser beam, overlapped to the 
trapping green beam, is focused by the same microscope 
objective far from the film surface. The observed portion 
of the red beam wavefront is approximately plane and the 
film reflectance than varies with cos[47rC(r)n/A] giving in- 
terference rings, n = 1.37 is the refractive index of the 
liquid mixture and A = 657 nm the red laser wavelength. 
The distorted film will then show ring shaped interfer- 
ence fringes centered on the trapped particle (Fig. 3b). 
The thickness profile around an isolated trapped parti- 
cle can be extracted (within an additive constant) and is 
reported here in Fig. 3a as open circles. A clean logarith- 
mic shape is found up to 200 /im, in perfect agreement to 
prediction in (|4]), as shown by the logarithmic fit curve 
(red line). 

To extract the force law one of the two traps is held 
fixed while the other is continously scanned through dif- 
ferent distances with a step of 2 /im. An image of the 
two particle is digitally recorded (Fig. 2b) for every 
scan step and subsequently used to extract particle po- 
sitions. Due to liquid drainage, the film is slowly but 
constantly thinning. When the thickness is so small that 
optical forces cannot balance the capillary interaction, 



one of the two particles jumps out of the trap and col- 
lapses onto the other. Until that time we can extract 
the attractive capillary force from the interparticle dis- 
tance. Calling k\ and k 2 the two trap elastic strengths, 
each particle will be displaced towards the other by a 
distance Axi = f/k{. The observed interparticle sepa- 
ration will be then smaller than trap separation by an 
amount: As = (Ax 1 + Ax 2 ) = f(l/h + l/k 2 ) = f/k'. 
Particle distance will actually fluctuate due to Brown- 
ian motion with a mean squared value given by 0: 
(Ad 2 ) = KsT/k' which can be used to experimentally 
determine k' = 29 pN//im. 




FIG. 4: (Color online) Intensity of the capillary interaction 
between two beads confined in a thin film, as a function of 
their relative distance. Open circles are experimental determi- 
nations. The black solid line represents a fit to the predicted 
power law (|12[) having the film thickness h as the only free pa- 
rameter. Red dashed line is the full Kralchewsky theory [12| 
for the same h parameter. The gray region marks the un- 
accessible range of separations s smaller than close contact 
distance 2a. Double log plot insert more immediately shows 
that the data are well represented by a power law. A 1/r law 
is also reported there, for reference, as a black dashed line. 

As already discussed, the logarithmic film shape in Fig. 
3 would lead, in the Nicolson approximation, to an at- 
tractive force decaying roughly as the inverse of interpar- 
ticle distance. More precisely we would expect to find 
the power law in ([12]) . Indeed a clean power law is found 
experimentally as evidenced by the double log plot re- 
ported in Fig. 4b. The overall behavior is very well 
fitted by the power law in ([T2|) leaving the film thick- 
ness at large distance h as the only free parameter. We 
find h = 2.2/im = 0.8a that makes an average gradient 
ip c = 0.12 confirming that the small gradient approxima- 
tion is justified for our system. For this fitted h value 
we get a power law exponent 1 — 2a = 0.86. The full 
Kralchewsky prediction for the same parameter values, 
also reported in Fig. 4b, is almost completely overlapped 
to our power law expansion, confirming the goodness of 
our approximations. 



We have shown how optical micromanipulation pro- 
vides an unprecedented tool for investigating capillary 
forces which govern aggregation and self-assembly of col- 
loids in liquid films. We provide a static force measure- 
ment of the capillary attraction between an isolated col- 
loidal pair and perform a direct test of the theoretically 
predicted power law. The exponent of the power law is 
found to be close to, but smaller than one, making cap- 
illary forces a quite unique example of very long ranged 
interactions in the mesoscopic world. The experiment 
opens the way to a variety of further developments ad- 
dressing the role of many body-effects, membrane elas- 
ticity, wetting properties, surfactant dynamics, hydrody- 
namic interactions in 2D. A deep insight into the na- 
ture of interface-mediated forces at the mesoscopic scale 
could suggest new routes to self-assembly of meso- and 
nano-structures [l9| . Optical trapping of colloidal parti- 
cles bound to lipid membranes [20[ could also provide new 
insights in the dynamics of biomembrane inclusions [21]. 
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